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Abstract. We establish Dahlberg's perturbation theorem for non-divergence 
form operators C = AV'^. If Co and Ci are two operators on a Lipschitz 
domain such that the Dirichlet problem for the operator Co is solvable for 
somep G (1, oo) and the coefficients of the two operators are sufficiently close in 
the sense of Carleson measure, then the LP Dirichlet problem for the operator 
Ci is solvable for the same p. This is an improvement of the Aoo version of 
this result proved by Rios in 1101 . As a consequence we also improve a result 
from [4] for the solvability of non-divergence form operators (Theorem l3.2ll 
by substantially weakening the condition required on the coefficients of the 
operator. The improved condition is exactly the same one as is required for 
divergence form operators C = div AV. 



1. Introduction 

This paper is a continuation of a long line of work, most recently advanced in 
[5 , on the solvability of the Dirichlet problem for elliptic operators with rough 
coefficients on Lipschitz domains with small constants. Here, we extend the results 
to non-divergence form operators satisfying a certain oscillatory Carleson condition 
which was left open in [4j. 

In particular, the non-divergence form results in 0] require a condition on the 
gradient of the coefficient matrix, since their results arise from considering the 
divergence form case first. They then change a non-divergence form operator to 
divergence form by allowing first order terms. 

Throughout this paper, the operators C which we consider are second-order, 
linear, uniformly elliptic and in non-divergence form. Precisely, C — a^^ {x)dij 
(we use here and throughout the paper the usual summation convention), where 
A{x) — (a'-' (x))^ . is a symmetric matrix with ellipticity constant < A < oo such 
that for aU x,^ e R", 

Me <eA{x)^<\-'\i\\ (1.1) 

We assume throughout that n > 2. 

The problem under consideration is the Dirichlet problem 

Cu = inD (1.2) 

u = g on dD, 
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where D C M" is a bounded Lipschitz domain. 

It is a fairly difficult task to even define the notion of a solution to the equation 
(jl.2p . Recall that in the divergence form case once can use Peron's method to 
construct a solution for coefficients that are merely bounded and measurable. This 
is not the case in our situation. For this reason we postulate what we mean by 
solving the continuous Dirichlet problem, denoted CD, following [10\ and [llj : 

Definition 1.1. (Continuous Dirichlet problem, CD) Given an operator C we say 
that the continuous Dirichlet problem is uniquely solvable in D (and we say CD 
holds for C) if for every continuous function g on dD there exists a unique solution 
u of (1121) such that u e C(D) n Wf;l{D) for some 1 < p < oo. 

From the results in |lj , if the coefficients a*^ are in VMO (see section [5] for the 
definition) and g G C{dD), then (11.21) has a unique solution Ug e C{D) n Wf^^{D) 
for all p, 1 < p < oo. By approximation, one can extend this result to allow 
for coefficients in BMO^o (also defined in section [2]) in a restricted range of p: 

i<p<po{qo) (c.f. m)- 

Definition 1.2. [L^ Dirichlet problem. Dp). We say that the Dirichlet problem is 
solvable for C in on D (or that Dp holds for £ on _D), for 1 < p < oo, if CV holds 
for £ and there is a constant C (depending only on C, A, n, D and p) such that for 
all g £ C{dD), the CD solution to (11.21) (which we shall denote by Ug) satisfies 

ll^%llLP(aD) < C\\9\\Lv{dD), (1-3) 

where N is the non-tangential maximal operator (see below). Unless explicitly 
stated, the assumed measure on dD is a, standard surface measure. 

Here the non-tangential maximal operator N (when necessary, N^) is defined as 

Nau{Q) = sup \u{x) \ , 

where ra{Q) denotes a truncated cone interior to D of aperture a based at Q G dD, 
i.e., 

TaiQ) = {xeD:\x-Q\<{l + a)Six)} n B,. (Q). 
Throughout, a > a*{D) > 0, with a* and r* > determined by the Lipschitz 
character of the domain and its size. Finally, 6{x) := dist{x,dD); Br{x) denotes 
the ball of radius r centered at x. 

In this paper we consider two fundamental questions. The first one is whether 
the solvability can be perturbed, that is, if Cq and Ci are two operators close in 
some sense, under what conditions does the solvability of the Dirichlet problem 
for one operator imply the same for the other? This question has a long history in 
the case of second order elliptic divergence form operators. For our purposes the 
papers f2j and [6] are of particular importance. Operators with first order terms 
are considered in [7]. 

The non-divergence form case has a considerably shorter history since new diffi- 
cult issues arise. In particular, the non-uniqueness of so-called weak solutions causes 
trouble in the most general case (see [9] and [13]). 

However, assuming, as we do, that the coefficients of the non-divergence form 
operators considered have a small BMO norm, one can establish the existence of 
strong solutions (i.e. solutions in Wf^^{D), c.f [1]). These are the solutions we 
consider in this paper. 
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The papers [TO] and [H] have made very good progress in setthng the question 
whether resuhs that hold in the divergence form case extend to non-divergence 
form operators. In particular, these papers show that if the elliptic measure of an 
operator Cq is in the Muckenhaupt Aoo{dcr) class, then so is the elliptic measure of 
an operator Ci under the same assumptions as in 6^. This implies that solv- 
ability of the operator Cq gives solvability of the operator Ci (for q potentially 
much larger thanp). The paper ^llj also considers first order terms (drift terms). 

In our Theorem 13.11 we settle the question of whether q can be taken to be 
the same as p, and the answer is affirmative if the coefficients of the considered 
operators are sufficiently close in the sense of Carleson measure. Analogous results 
for divergence form operators have been established before ([2] and [6]). We do 
not consider first order terms as JTl] does; however, this can be done, as all the 
necessary ingredients are already in place. We choose not to do it here to make our 
already very technical exposition more readable. 

The second fundamental question we settle here is the question of finding a broad 
condition on coefficients of the non-divergence form operator that guarantee 
solvability. Again, the case of divergence form operators serves as a model. There 
are two particulary important results to mention here. In [8] it was established 
that under the assumption that t|VA| is a Carleson measure, the LP Dirichlet 
problem is solvable. In 4, this condition was relaxed (the gradient is replaced by 
an oscillation- type condition), and it was also shown there that given p £ (1, cxd) the 
LP solvability depends on the norm of the Carleson measure. If the norm is small, 
the LP solvability for particular p holds. Moreover, since first order terms are also 
considered, the results in [8 and 4^ do apply to non-divergence form operators, but 
only under the stronger gradient condition t \^A\ . 

The missing piece to prove this under the weaker "oscillation condition" is a 
strong perturbation theorem for non-divergence form operators which we establish 
here. Hence the result in [1] is substantially improved. This is formulated in 
Theorem 13.21 As previously mentioned, the weaker A^o version of this result is 
already done in [Tl] . 

The structure of this paper is as follows: section 2 contains a few definitions 
needed to formulate two main results, which is done in section 3. Section 4 expands 
the list of definitions and introduces a few technical preliminaries. Section 5 contains 
the proof the the perturbation result and section 6 is dedicated to the LP solvability 
under the Carleson condition on the coefficients of our operator. Finally, sections 
7 and 8 contain the proofs of two auxiliary lemmas. 



where Je is the average value of / on E. Then / G BMO(M") (i.e., / has bounded 
mean oscillation) if 77 G L°°(M+,M"). Moreover, ||/||bmo ~ II'7/I1l°=(r+,R'»)- 

Let 7y(r) = ||ry(r, •)||ioc(R„). We say / G VMO(R") (/ has vanishing mean os- 
cillation) if limr^o+ vir) = 0. Finally, a function / G BMOg(R") if r]{r) G ^{q), 
where $(^?) is the collection of all non-decreasing functions 77 : M"*" — ^ M"*" such that 



2. Basic definitions 



Given / G Ll^{R"), let 
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there exists a C > such that ri{r) < g for ah r < C,. These spaces (BMO, VMO, 
BMOg) can be restricted to a Borel set G using standard methods. 

The setting for our work is a Lipschitz domain D. A bounded, connected domain 
D e M" is called a Lipschitz domain if there is a finite collection of right 

circular cylinders and Lipschitz functions 0^ ((/>, : M"^^ — > M, and there is an 
L > such that for all x,y E M"^-'^, \(j){x) — (p{y)\ < L\x — y\) such that the 
following hold: 

(i) The collection of cylinders {li} covers the boundary, dD, of D. 

(ii) The bases of the cylinders have positive distance from dD. 

(iii) Corresponding to each pair {Ii,(f>i), there is a coordinate system {x,s) with 
X e M"^^, s e M such that the x-axis is parallel to the axis of Ii and such 
that n Z? = {(x, s) : s > (l>i{x)} H h and hCidD = {{x, s) : s = (t)iix)} f] h. 

Without loss of generality, we will assume that D is containted within the unit 
ball centered at the origin of M" and that D contains the origin, i.e., we assume 
D C Bi{0) with e D. 

For points Q E dD and r > 0, we denote the boundary ball of radius r at Q 
by Ar{Q) = Br{Q) n dD. The Carleson region Tr{Q) above /S.r{Q) is given by 
Tr{Q) = Br{Q) n D. We say that a measure /i on D is a Carleson measure if there 
is an Af < cxd such that 

sup , , = M. 

r>0,QedD Cr[Ar[Q)) 

3. Main results 

The following perturbation theorem is modelled on Dahlberg's theorem (Theo- 
rem 1 in [2], re-proven as Theorem 2.18 in 6 ). 

Theorem 3.1. Consider operators Cq, Ci, with Ck = a^^{x)dij on a Lipschitz 
domain D, e{x) = laQ{x) — a^-'(2;) j and sl{x) — sup^^^ ^^^^ (j,) |£(2)|. Let A > 
be the ellipticity constant of the operator Cq, and let 

/ TTV dx^eo< oo. (3.4) 

Assume that the L^ Dirichlet problem is solvable for the operator Cq with a constant 
Cp > in the estimate (II. 3p for some 1 < p < oo. 

There exist constants go > (independent of p) and M = M{p, D, X,Cp, g) > 
such that if £ BMOg with g < go, and if Eq < M , then the L^ Dirichlet problem 
is solvable for the operator C\. 

Remark 1. This theorem is a direct improvement of Theorem 1.1 in 10 and The- 
orem 2.1 in [TT] where a statement of the type "Dp for Cq =^ Vg for Ci was 
established with q » p. 



Remark 2. It suffices to assume that the condition Eq < M in Theorem 13.11 onlv 
holds for all Carleson regions such that r < rp for some > 0. This is due to 
the comparability of the elliptic measures of two operators whose coefficients are 
the same near the boundary (see Lemma 2.15 of [TO]). 

Remark 3. The theorem can be formulated on more general domains. In fact, we 
never explicitly use that the boundary of D has a graph-like structure. The minimal 
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geometric structure needed is that 13 be a chord-arc domain and non-tangetially 
accessible. 

Remark 4. The number qq is chosen such that if a'-' G BMOg for any g < qq, CD 
holds for the operator C = a^^dij on D and so that the value of po(ffo) referenced 
in the discussion after Definition 1 1 . 1 1 satisfies po(f?o) > 2. 

Theorem 3.2. Let I < p < oo, let < X < oo be a fixed ellipticity constant, and 
let D be a Lipschitz domain with Lipschitz constant L. Let C = a^^dij be an elliptic 
operator with ellipticity constant X. 



is the density of a C'arleson measure in D with Carleson constant M , then there 
is a constant C(p, A) > such that if L < C'{p,X) and M < C{p,X), the Dirichlet 
problem Dp is solvable for the operator C. [Here avg(a*^(z)) is the average of the 
coefficient a*-' over the ball Bg(^z)/2{x)-] 

This is a substantial improvement of Corollary 2.5 in 4 , in the spirit of Corollary 

2.3 in the same paper for divergence form operators. This also improves Theorem 

2.4 in TT. Note that we do not consider drift terms as lHj does. 



To enhance the readability of this paper, we have kept our notation the same 
as in [10 . We rely heavily on certain results from [10 in the technical part of this 
paper. 

Throughout the paper, we use A < _B to mean there is a constant C, depending 
on, at most, n, A,?] and D such that A < CB; similarly ioi A> B. li A < B and 
A> B, then we say ^ « B. 

Let C be an elliptic operator for which CD holds. By the maximum principle, 
the mapping g i— > Ug{x) is a positive linear functional on C{dD) for each fixed 
X € D. The Riesz representation theorem then gives a unique regular positive 
Borel measure on dD such that 



This measure is called the harmonic measure for C on D. 

Given a non-decreasing function rj, we denote by 0{X, rj) the class of operators 
£ — a^^dij with symmetric coefficients satisfying (jl.ip such that a'-' G BMO(M") 
with BMO modulus of continuity rj in D. We use 0(A) if there is no restriction on 
the regularity of £. 

The theory of weights plays an important role in what follows. Given a p, 
1 < p < oo, and two measures fi and v on dD, if is absolutely continuous with 
respect to i^, let fc = Then, we say that /j, G Ap(dv) if there is a constant C < cxd 
such that for all boundary baUs A (i.e., for some r > 0, Q G dD, A = Ar{Q)), 




4. Notation and Technical Preliminaries 
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We say that G RHp{dv), the reverse- Holder class, if there is a constant C such 
that for all boundary balls A, 

^ k^dv] <C-j—- I kdv. (4.6) 



A 

/ — P 



Note that /i £ Ap{dv) if and only \i v E RHpi{dii), with p' = The best 

constant C in (14. 5p is called the Ap "norm" of fj, and is denoted Ap{^\di^). Recall 
that the assumed measure on dD is cr, standard surface measure, so by /i € Ap, we 
mean /x S Ap((iCT). Also, these classes of measures (or weights) are related: 

y RHp,{dv) = y Ap(d;/) Aoo(d;^). 
p'>i p>i 

A crucial ingredient in what follows is the fundamental theorem relating weights 
to solutions of elliptic partial differential equations (first proved by Dahlberg in [2] ) : 

Theorem 4.1. Let uj be the harmonic measure with respect to C on D , and let /i 
he a Borel measure on dD. Then the following are equivalent: 

(i) Lu e Aooidfi). 

(a) There is a I < p < oo such that Vp holds, i.e, 

ll^"gllLP(aZ?,(iM) ^ Cp\\g\\LP(^gD^dfj,}- 

(Hi) UJ is absolutely continuous with respect to fi and uj G RHpi{dfj,) (where, again. 

Lemma 4.2. (Theorem 2.5 in |10j ) // we let p G [n, oo), w G Ap, then there is a 
constant Qp such that if r] € *&(pp) and C G 0{X,r]), for any f G L^{D,w), there 
exists a unique u G C{D) D W^'P{D, w) such that Cu — f in D and u — on dD. 



Then, with / and u as in Lemma 14.21 for each x & D, the maximum princi- 
ple implies that the positive linear functional / — )■ —u{x) is bounded on U'{D). 
The Riesz representation theorem then gives us the unique non-negative function 
G{x, ■) G LP {D) {p' = ^) such that 

{.x)^- I G{x,y)f{y)dy. 



JD 

This is the Green's function for C in D. 



Definition 4.1. Given C G 0(A), v G L\^^{D) is an adjoint solution of Cva D if 

vC(j) dx = 



for aU (f) G C^{D). In this case, we write C*v = 0. 

Recall that we are assuming D C i?i(0). Pick a point x G dBg{0), and let 

Piy) = G£.Sio(o)(S,2/) in Bio(O), (4.7) 

where Gc,Bio(.o) is the Green's function for C in _Bio(0). The following technical 
estimate is quite useful. 
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Lemma 4.3. (Lemma 2 in 5 ). Let G{x,y) be the Green's function in D for 
C G 0(A). Then there is a constant ro depending on the Lipschitz character of 
D, such that for all Q G dD, r < ro, y € dBr{Q) Pi Vi{Q) and x ^ T^riQ), the 
following holds: 

with p as defined in ^Ji), B{y) Bs[y)/2{y) and p{B{y)) = /g^^^ p{y)dy. 

Following Rios, '101, we define a modified area function and non-tangential max- 
imal function which are adapted for the non-divergence form situation. 

Definition 4.2. [Area functions] For a function u defined on D, the area function 
of aperture a, SaU, and the second area function of aperture a, AaU, are defined 
as 

Sau{QY ^ [ frfrlx \Vu{x)\'^ p{x)dx, and 

Jr^(Q) P{B{x)) 

A^u{Qf^ [ firlx \y^u{x)\^ p{x)dx, 

with p{x) as in (14. 7p . B{x) = i?5(x)/2(2;), and Q G dD. 

We also recall Rios' modified non-tangential maximal function, 

(n^{v)) {Q):= sup / v{yf^^dy. (4.8) 

Here, Bo{x) denotes B{x, ^^). 

5. Proof of the perturbation theorem 13.11 



The structure of our proof owes much to the proof of Theorem 2.18 in [5]. We 
use Sa, Aa and as defined in the section above. 

Let Co and Ci be two operators as in Theorem 13. 1[ and consider any con- 
tinuous boundary data g. We first establish that CV holds for Ci. We observe 
that \\Ao - AiWl^f^^-^ ^^o- Since G BMO^ and holds with £o < M, 

aY G BMOg+£, where e can be arbitrary small (it depends on M). So if M is 
made small enough, we can ensure that g + e < qq. Recall that go is chosen to 
guarantee that CD holds for any operator Ck = a^dij, as long as a^^ G BMO^ and 
g < go- 
Notice also that if A is the ellipticity constant of Co, one can guarantee that the 
ellipticity constant of Ci stays bounded away from zero, say by A/2, by making M 
smaller if necessary. 

Hence we can talk about solutions uo and ui to the corresponding Dirichlet prob- 
lem with the same boundary data g for Co and Ci, respectively. Let F — uo — ui. 
If follows that CoF = —Coui, so 



F{x) / Go{x,y)CoUidy. 



D 



Here Go is the Green's function of the operator Cq. 

We will use the following two lemmas and defer their proofs until later. The 
following is analogous to Lemma 2.9 in [6]. 
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Lemma 5.1. There exists a constant C = C{X,n) such that under the hypotheses 
of Theorem \3.1\ 

NF{Q)+N{5\VF\){Q) < CeoM^,{A&^i){Q)- 

The next lemma is analogous to Lemma 2.16 in [6]. 

Lemma 5.2. Let a > 0. Then there exists < /3 < a depending only on the 
dimension, the number a and the Lipschitz constant of the domain D such that the 
following holds: 

Suppose that Si3{F)(P) < A for some P in a surface ball A = A{PQ,r) C dD. 
Then there exists c > 0, S > depending only on the Lipschitz character of the 
domain D and the ellipticity constant of the operator Cq such that for any 7 > 

a({QeA ; ^^(i^^) > 2A, iV„(F) < 7A, 7V„((5 1 VF|) < 7A, 

iV„(F)A,(ui) < (7A)2}) < C7V(A). (5.9) 

Assuming Lemmas 15.11 and 15. 2[ we have 



N{F)Pda< f (^N{F)P + N{6\VF\)P'^ da 



dBi JdBi 



<Ceo / {ALM»ui))Pda 

JdBi 

<Ceo [ {M^SA^ui))Pko^dcuo 

JdBi 

<C'eo f A^{ui)Pk^Uujo, 



where fco = Recall that fco G RHpi{da) (since we are assuming Vp solvability 
for Cq) is equivalent to /cq"^ G Ap(rfwo). 
We then have 

/ (n{F)p + N{5\VF\)P) da < Ceo [ Aaiui)P da 

JdBi ^ ' JdBi 

< C'eo / Sc5(ui)P da 

JdBi 

< C'eo f Sfs{ui)Pda (5.10) 

JdBi 



< C'eo / {Sp{uo)P + Sp{F)P)da, 

JdBi 

where we used Theorem 2.19 from [10] and also Theorem 2.17 of [10] for 
II 'S'cq(ui) II LP ~ ||S'/3(ui)||lp. It remains to deal with these terms. 

First, we note that /^^^ {Scdi{uo)Y da < Jgg_^ f^da, using Theorem 2.17 from 
[10) and our assumption that the Dirichlet problem is solvable for Co- 

For the second term, we will use the good-lambda inequality from Lemma 15.21 
more specifically its Corollarv l8.ll According to it, we have an estimate 



S0{F)Pda<2C (NiF)P + N{6\yF\)p)da+ S0{uo)p da. (5.11) 

IdBi JdBi ^ ' JdBi 

The term Jgg_^ S'^(uo)^ da is harmless and can be estimated by Jgg_^ fP da, as above. 
Now we put (|5.10p and (|5.1ip together to obtain 
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/ (n{F)p + N{S \VF\y) da < Cso I (n{FY + N{5 \VF\Y + f) da. 
Hence for Eq sufficiently small so that Ceo < 1/2, we have that 
/ {n{FY -VN[8\VF\y\ da <C \ P da. 

JdBi ^ ' JdBi 

This is the estimate required for the solvability of the Dirichlet problem Vp for 
£i, as ui = uq—F, and for uq we have the needed estimates due to the solvability 
for Cq. Since for a' < a we have a pointwise estimate: 

Na^'UiiQ) < N^uiiQ) + N^{5 \Vu^\){Q), 

the theorem follows. 

6. Proof of Theorem 13.21 

Now that we have established the perturbation theorem, we can easily dispense 
with the proof of Theorem [321 The first part of the proof, dealing with the smooth 
perturbation of A, is exactly the same as the smooth perturbation part of the proof 
of Corollary 2.3 in [4]. We repeat it here for convenience. 

We prove this in the flat case; the general result will follow from a change of 
variables by Necas and Stein (see, e.g., p. 2 of [4i for details). The notation avg(a) 
at a point (j/, s) represents the average of a over the ball of radius s/2 centered at 
{y,s) (denoted Bg/2{y, s)). Given a matrix coefficient a{x,t) in R", set d{x,t) — 
J a(u, s)(j)t{x — u, s — t) ds du, where is a smooth bump function supported in the 
ball of radius 1/2 and (ptiu, s) — t'^-ipd/ft, s/t). 

We are assuming that 

(sup{|a(2/,s)-avg(a(a:,0)|' : {y,s)eBt/2{x,t)Y) ^ (6.12) 

is a Carleson measure with small norm. 
We aim to establish three facts: 

t\Va{x,t)\^ dxdt (6.13) 

is a Carleson measure with small norm, 

(sup{|a(y,s) - a(j/,s)|2 : (y, s) e B./^ix^t)]) ^ (6.14) 

satisfies the hypotheses of Theorem 13. 1[ and 

a{x, t) e BMOg for g = q{M), with ^ as M ^ 0. (6.15) 

Given the results in U , the condition (|6.13p implies that Dp holds for the operator 
with coefficients A. Using (|6.15l) . if M is chosen sufficiently small we will have 
g < go- Combining this with (|6.14p . as in the proof of Theorem 13.11 above, yields 
that A is in BMO^, for some g < g^. Thus, the hypotheses for Theorem 13.11 are 
satisfied and, therefore, Vp holds for the operator with coefficients A. 

That (16.131) follows from the hypotheses of Theorem 13.21 is a straightforward 
calculation; apply the gradient to (j)t{y,s) and subtract a constant from the Uij 
inside the integrand to see that 

\Vdix,t)\ < Ct-^ (sup{la(y,s) - avg(a(x,i))l : (y, s) £ Bt/2{x,t)]) . 
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The proof of (|6.14p is equally straightforward; add and subtract the constant 
avg(a(a;, t)) inside the difference. For precise details see [4] and a similar calculation 
in [10]. 

It only remains to prove (|6.15p . Choose an arbitrary point {x,t) in our domain. 
We shall check that the function a is BMO near this point. Consider a ball B of 
radius s > centered at the point {x,t). There are three cases to consider: 

(i) a small ball, with s <t/2, 
(u) a large ball, with s > 2t, 
(iii) an intermediate ball, with t/2 < s < 2t. 

As we shall see, only the cases (i) and (ii) are fundamental. Case (iii) is merely 
a combination of the approaches taken in (i) and (ii). 
In case (i), (|6.12p trivially gives that 



sup |a(j/,u) - avg(a(a;,t))| < M, 



hence 

oscBj/g(x.t)a = max 



sup a^^{y,u)— inf a^^{y,u) 

(y,u)&Bt/s(x,t) (y,«)eBt/8(a;,t) 



< Ml/2, 



From this 

oscB^(^.t)a < oscB,,2{x,t)a < M^/^ (6.16) 

as the ball or radius t/2 can be covered by a fixed number (depending only on 
dimension) of balls of radius t/8. This immediately gives 



\Bs{x,t)\- 



Bs{x,t} 



\a'^{y,u) - avgsJa^J')! dy du < Af^/^ 



hence g < CM^'^. 

If (ii) holds then Bs{x,t) intersects the boundary {t = 0} at a large set of area 
of order s"^^. One might think of 13 n Bs{x, t) as a subset of a larger Carleson box 
T(A), where A is a surface ball on the boundary {t — 0} of radius comparable to 
s (a multiple of s where the constant depends on the dimension of our domain). 
Therefore, it will suffice to prove that 



T(A) 



a'->{x,t) - avg7,(^)(a*-') 



dxdt<M^/^s'', 



from which again g < CM^/"^ on such balls. 

In fact, the exact average that gets subtracted off in the BMO norm does not 
matter, so we might as well prove that 



/T(A) 

Here Tfc(A), k = 0,1,2, 



a'\x,t) -avgy^(^)(a*J) 



dxdt<M^''^s''. 



is defined diadically by 



Tfe(A) = {(a;,<)eT(A);<e(2- 



-fc-i. 



It follows that T(A) is a disjoint union of rfe(A), fc > 0. 
By (|6.16p we immediately get that 

a''{x, t) - avgj,„(A)(a'J') < M^'\ for aU (x, t) e To(A). 
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Now consider {x,t) G Ti{A). By using (|6.16|) twice we get that 

a'^ix,t) - avgj,^(^)(a^^)| < 2M^/^, for all {x,t) e ri(A), 
and inductively 

a'Hx, t) - avgj,„(^)(a'^') < {k + l)M^'^, for all {x, t) G Tk{A). 



Hence, 



T(A) 



a'\x,t) -avg-r^(A)(a'^) 



dx dt 



a''{x,t) - avgy^(^)(a'^) 
< ^(fc + l)Af i/2(2-'=s") « Afi/^s", 



since |Tfe(A)| « 2-'=s". 

Case (iii) is a combination of these two approaches where one considers the 
integrals on pieces Bs{x,t) n {{y,u)- u G (2^'°^^i, We leave the details for 

the reader. 

By combining (i), (ii) and (iii) we see that G BMOg for q < CM^/^, where 
M is the bound on the Carleson measure of (I6.12p and C > is a constant that 
depends on the dimension of our domain D. 



7. Proof of Lemma EH] 
Proof of Lemma \5.1[ From [lOj, Lemma 3.2, we know that 

Nc.F{Q) < soA'L„iAs,ui)iQ) (7.17) 
for some a slightly larger than a. We will also show that 

(7V„((5|VF|))'(g) (7.18) 

< Ns,{F)iQ)Na{S |VF|)(Q) + (iV5(F))' (Q) + eoNaiF){Q)Aa{ui)iQ). 

Combining these two yields the lemma. Thus it remains to show (j7.18l) . 

To this end, we fix Q G dD, x G TaiQ). Also, find the required value for vq 
in Lemma 2.14 of [10] (if necessary, making < j), and then choose r* < ro/2, 
where r* is the truncation level of F^iQ)- 

Under these assumptions, if we take y e = B{x, ^^) then y G Ta{Q) for a 
slightly larger cone and also y G dBr{Q) for r < rQ. Hence, Lemma 2.14 of [10] 
can be applied to all of the points in our integral. Lemma 2.14 in |10] provides the 
estimate: 

Go{0,y) ^ 6Hy)p{y) 
wo(A.(g)) p{Biy)) ' 

where r = \y — Q\. For y G Bq, we have that ^S{x) < r < ^S{x). We observe 
that for r in this range, all values of ijJo{Ar{Q)) are comparable to the value of 
a;o(A^(a,)((5)), as the measure is doubling. Also, let 6 := S{x). 
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Following [6], we start with 



where S,- = -B(a;, r). However, we need to average twice. Thus, we estimate 

(7.19) 

with Q^{s) = {Pi - \)s + f , g2{s) = (/?2 - \)s + f , with /3i < i < i < /?2. The A's 
are yet to be determined. 
Then, 



< 



Here Ao is the matrix of coefficients {a^)- We first estimate the contribution to 
()7.19p by Ii. Integration by parts twice yields: 



IiLOoiAsiQ)) ^ CoiF')GoiO,y)dy 

= V / (Goa}^d,{F^)i^, - d,{al'Go)F^>^,) da + [ C*{Go)F^ dy. 

JdB^ ^ ' Jb^ 

However, Cq{Gq) — 0, so we are only left with the two boundary terms. Hence, 



22 (s) 

hdr 



I y 



^ [Goal^d,iF'):^,-d.{al^Go)F'iy,)dy 



Goad' {dj{F^)u, + d,{F^y, + F^d^u-j) dy 



/ G qQq F^Vi Uj da 
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We estimate this term by term, starting with 
1 



SujoiAsiQ)) 
1 



E 



-Se2(»)\-^ei(s) 



Goal' {d,{F'y, + d,{F^y,] 



dy 



< 



|i.||VF|i^ 



dy 



It 



And 



It ds < 



F\\^F\^M.dy 



<N^{F)Ns,{5\VF\). 



(7.20) 



Recall that a must be chosen a little larger than a. The parameters Pi determine 
the size of a, as we want all points in Bp^j C Ta{Q)- The this choice is irrelevant 
as long as Ta{Q) is still a nontagential cone. 

Next, we look at z/, the outward unit normal for Br arising from our first inte- 
gration by parts. We know Vj — ^ when |x| = r, and di{vj) = '^j^- Thus, for 

X G Bq2{s)\Bqi(s)i Pi^ — \x\<, /32S whence \diiyj\ < |. This leads to 



SujaiAsiQ)) 



E 



< 



GoaiF'd,{v,) 



dy 



1 



F' 



-Be2(s)\-^8i(s) 



p{B{y)) 



dy := ll 



and 



1 



Il< 







F^^^dy<(N,iF)^; 



p{B{y)) 



For the last term, we see 
1 



E 



< 



1 



< 



r^2 piy) 



da ds 
GoF^dy 



{B^^s\Bs/5)U(Bs/6\Bi3-^s) 

■ dy < (^Ns.{F)y . 



.... piBiy)) 

We now turn to estimating by the fact that uj is doubling. 



\h\ < 



/ 1^1 l^of I ^^ff^,, dy<eo[ |F||VV 



5^{y)p{v) 
PiBiy)) 



(7.21) 



(7.22) 



dy. (7.23) 



Here we are using the fact that CqF — —CqUi — — (£i + e^^ dij)ui = —e^^dijUi and 



sup^es, \£'Kz)\ < sup^es, 
IS-H we get that a(a;) < £o- 



|e(z)| — a(x). Using the condition from Theorem 
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Thus, 

1 



/ Ihl drds<- I/2I dr (7.24) 

Jeiis) " JpiS 



<eo f \F\\y^u,\^^^^dy<eoN^{F)AM 
Jb,,, P(B{y)) 



By combining ((71^ . ([7:^ . ([7211), ([7:^ and ((TTMl) . the lemma is proven. □ 

8. Proof of Lemma [Q] 

We now prove Lemma 15.21 the good-lambda inequality which is crucial for esti- 
mating S{F). 

Proof of Lemma \5.'d Let us call E the set 

S = {q e A; SpiF) > 2A, N^{F) < 7A, 

N^{S\\7F\)<-fX,N^iF)AM < (7A)'}. 

It is sufficient to prove that uj{E) < Cj'^uj{A), since we already know that oj £ 
Aao{d(j). Standard arguments (see [3] or [5]) show that since S 13(F) > 2 A on E, 
we can choose 7 > sufficiently small so that Sp^rriF) > A/2, where each cone is 
truncated at height rr for some fixed < r < I. It follows that 

co{E) <^J^ Sl,,{F){Q)doj{Q). (8.25) 

We will introduce the following notation. Let 

2^ - U ^P-.rriQ)- (8-26) 

For a' G {I3,a), let us consider a smoothed-out version of the set [Jq^e^oi' ,Tr{Q)- 
We denote by Va' the set with the properties: 

(ii) dT>a' is smooth except at E and \Vi'{Q)\ < C /5o{Q) for Q G dVa', 

(iii) Pa' C Pa" if a' < a". 

Here is the outer normal at the boundary and 6 — So denotes the distance to the 
boundary of the original domain D. We now work with (|8.25p . 

^{E) < ^f(f \VF\'^^^dx]d^{Q) 

^ ^fif NFfGo^^] dcoiQ) (8.27) 

< -S / \VF\''Godx 



< {AoVF-VF)Godx. 



Here ^ {Q e dD; x G r^(Q)} and a' G (/3, a). Now, 

AoVF ■ VF = Co{F'^) - 2i^i;oi^, 
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SO there are two terms to estimate 

{AoVF ■ VF)Go dx 



= [ CoiF^)Godx- [ 2FCoFGodx. 



(8.28) 



Let us denote these two terms by I\ and l2- We first deal with I^. Recall that 
£0^0 = —^(0), hence integration by parts gives us only two boundary terms 



h < f 4^ diiF^)GoUj da - [ dj{c^^Go)F'^Uida. (8, 



29) 



Note that, strictly speaking, these two boundary terms are not well-defined. To fix 
this, we again use the averaging technique introduced before. We integrate over 
the interval [a', a"] C a) and get instead solid integrals 



h < c 



+C 



c^^di{F'^)Go5-'^Vjdx 

dj{ai^Go)F'^5-'^Uidx 



Now, for simplicity let V — Va" \ V^' ■ We see that 



Ii < C ( \F\ \VF\GQ5-'^dx 

" I il 
- / (/ 

I il 



(8.30) 



dx 



JF||Vi^|Go- 



duj{Q) 



jF\\VF\5j^dx\du^{Q). 



< 



1/2 



■ ' piB{x)) 



dx 



X 



(8.31) 



'1 



rf,iQ)nv 



piBix)) 



\ 1/2 

dx ) du{Q). 



The key is that if a; e r^(Q) n V then x € TdQ') for some Q' € E and the set 
r^((3) n 2? is of diameter proportional to 6{x) and its distance to dD is also of 
6{x) size. This implies we can control the two solid integrals the the last line by 

N„{F){Q')N„{S\VF\)iQ'). This gives 



since the measure lo is doubling. To estimate l\ we integrate by parts one more 
time. We get 

^'""'■^ ' [ ai^GoF^iyiiyj5-^da 

J&D 



It < C 



a^Godj 



V 



dx 



+ c 



(8.32) 
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The first term of (|8.32p will give us two additional terms, depending on where the 
derivative dj falls. By the chain rule, 



< 



CIFIIVFI CF2 



52 



Here we use the fact that the real distance function 5 can be replaced by a smooth 



distance function so that \\7S 



S 2 and also iVvil < C6 ^. Hence, the first 



term is of the same type as , and the second one can be bounded by 



\Ff ^J^^^dx] dLj{Q). 



p{B{x)) 



5.33) 



da < CYX^uj{A) 



Thus, this term can be dominated by C Jq^2A (^■^aiF){Q') 
Finally, (|8.32l) has one additional boundary term, which again has to be averaged 
out. So we need to use the wiggling technique one more time. Without going into 
too much detail, this will again turn the surface integral into a solid integral over 
a set we call V (essentially of the same type as D): 

F^iy,iy, 



Oq Go - 



-dx 



< 



F^Gn6-^ dx. 



Notice that this term is similar to (I8.33p . so the same estimates can be applied. 
This establishes 

i/ri < c-f^x'^LJiA). 

Now we deal with l2- As before, we use CqF — —e'^^dijUi, where e^^ — — . 
This gives 



< C 



I e{x) \F\ jV^uij Gadx, 



.34) 



where e{x) = max|£y (a::)|. We turn this back (by Fubini) to into two integrals 



< 



e{x) \F\ IV^uil 



Qe2A \''T^,{Q)r\V^ 



Go 
w(Aj, 



■ dx da{Q) (8.35) 



£{x) \F\ \V^u 



Qe2A \ "'r„,(Q)nX'„ 



' PiBix)) ) 



da{Q). 



By Holder: 



h < 



s{xr\F\' 



1/2 



Qe2A \Jr^,{Q)n-D^ 



p{B{x)) 



■ dx 



5.36) 



1/2 



r„,(Q)n-D„ 



p{B{x)) 



dx 



da{Q). 



As a' < a, it can be arranged that either Ta'iQ) n Va' C TaiQ') for some Q' ^ E 
or N4F){Q)Aa,{ui){_Q) < i-f\)\ 

Indeed, if Q € then the fact that iVa(F)((5„)Aa(wi)(Q„) < {'jX)'^ for a 
sequence of Qn £ E converging to Q implies the same for Q. In this case we just 
take Q' ^Q. 
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Otherwise d = dist{Q,E) > 0, and this gives that Ta'iQ) fl Va' only contains 
points of distance S > d. Hence by making a sufficiently large we will have Ta' {Q) n 
Va' C Ta{Q') for all points Q' ^ E such that dist(Q, Q') w d. If follows that 

\Jr^,{Q)nv^, p{B{x)) J 

On the other hand, 

)l/2 
<Esrr{Q)Na{F){Q'), (8.37) 

where 

/ 



(suPB(x,5(a;)/6)e(a;)) 

AQ) 




To see (|8.37p we cover the set To,' (Q) H Va' by a union of balls of diadic diameters 
fc e Z, with each such ball of approximate distance to the boundary such 
that each point x G Ta'{Q) H Va' belongs to at most K balls. (Simple geometric 
considerations imply that K will only depend on the dimension, the number a' and 
the Lipschitz constant of D). On each such ball, the square of the solid integral 
on the left-hand side of (|8.37|) can be estimated by C (sup^^gBi ^{x)'^) Na{F)^{Q'). 
After we sum over all the balls we get the expression CKE^^^{Q)Na{F)'^{Q'). It 
follows that 

h < C I AaM{Q')Na{F){Q')E3rriQ)dwiQ) 
JQe2A 

< C-r^X^ [ EsrriQ) dujiQ) (8.38) 

JQG2A 

1/2 



\jQe2A 



El^{Q)duj{Q) < Cio{2A) 

Qe2A 

by Rios's work (see p. 683 of [10]). Note that this is the only place we are using 
p.4p . and we do not use the fact that it is small, only that it is finite. This 
establishes the good-lambda lemma. □ 



Corollary 8.1. Lemma \5.^ implies that for any 1 < p < oo: 

[ S{FYda<C{q)( {{N{FY + N{5\VF\Y)da+ ( Siuofda, (8.39) 

JdD JdD ^ ' JdD 

where the square function S is defined over cones of smaller aperture than the 
modified nontangential maximal function N. 
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Proof. Indeed, the Whitney decomposition and Lemma 15.21 gives us 



dD 



SfiiFfda < C 



dD 



NaiF)P + N^{6\\7F\)P]da 



+ / iA^{u,)N^iF))P/^ da 

JdD 



for some (3 < a. This imphes that for any e > 0, 

[ Sf}{FYda <C{e) ( (No,{F)p + NaiS \VF\)p) da + e [ Ao,{ui)p da. 

dD JdD ^ ' JdD 

By Theorem 2.19 of [TO], since ui is a solution to CiUi = we have pointwise 
estimate Aa{ui) < CScaiui) for some c > 1 depending only on the dimension n. 
Also by [S] (see also Theorem 2.17 of [TO]) for solutions we have ||5cq(ui)||i,p < 
C||iS'^(mi)||lp with C only depending on the ellipticity constant, the numbers ca 
and /3 and the dimension. 
This gives 



dD 



'fi{F)Pda <C{e) [ {n^{F)p + N^{5\VF\)p)da + Cie [ Sp{ui)P da. 

JdD ^ ' JdD 



We can write Sfi{u{)P < C2{Sp{uo)P + Sp{F)P). Choose e so that CiCae < 1/2 
(this allows the term CiC2sSi3{F)p to be incorporated into the right-hand side). It 
follows that 



dD 



Sfs{F)Pda < 2C{e) 



dD 



N^{F)P + N^iS\VF\)P da 



dD 



Sp{uo)P da. 



□ 
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